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It is shown that the statistical properties of fully developed 
hydrodynamic turbulence are characterized by an infinite set 
of independent anomalous exponents which describes the scal- 
ing behavior of hydrodynamic fields constructed from the sec- 
ond and larger powers of the velocity derivatives. The energy 
dissipation field e(t, r) and the square of the vorticity are the 
simplest examples of such fields. A physical mechanism re- 
sponsible for anomalous scaling is discovered and investigated. 
We call this mechanism the telescopic multi-step eddy interac- 
tion. The essence of this mechanism is the existence of a very 
large number (R/r]) Aj ^ 1 of channels of interaction of large 
eddies of scale R in the inertial interval with eddies of viscous 
scale n via a set of eddies of all intermediate scales between 
R and r\. The description of this mechanism in the consistent 
analytical theory of turbulence based on the Navier Stokes 
equation in the quasi Lagrangian representation is presented. 
In the diagrammatic expansion of the correlation function of 
the energy dissipation field K ££ (R), we have found an infinite 
series of logarithmically diverging diagrams. Their summa- 
tion leads to a renormalization of the normal Kolmogorov-41 
dimensions. For a description of the scaling of various hy- 
drodynamic fields an infinite set of primary fields O n with 
independent scaling exponents A n was introduced. We have 
proposed a symmetry classification of the fields O n enabling 
one to predict relations between scaling the behavior of dif- 
ferent correlation functions. For instance the principal con- 
tributions to the irreducible correlation functions of all scalar 
fields constructed from velocity derivatives possess the same 
scaling behavior with a so-called "intermittency exponent" \i. 
Further we formulate restrictions imposed on the structure of 
correlation functions due to the incompressibility condition, 
e.g., the simultaneous correlation function (ev) (where e is 
the energy dissipation rate) is equal to zero. Experimental 
test for the conformal symmetry of the turbulent correlation 
functions are proposed. It is demonstrated that the anoma- 
lous scaling behavior should be revealed in the asymptotic 
behavior of correlations function of velocity differences. A 
way to obtain the anomalous exponents from experiments is 
described. 



INTRODUCTION. 

Problems related to statistical properties of fully devel- 
oped hydrodynamic turbulence have been intensively at- 



tacked during the last century. Nevertheless no consensus 
has yet been achieved on the behavior of the correlation 
functions of the turbulent velocity. Investigators agree 
that they should be treated in terms of a scaling behav- 
ior but the character of this scaling is still under discus- 
sion. An analytical theory of scaling behavior of hydro- 
dynamic fields constructed from the velocity derivatives 
is presented in this paper. We begin with a brief overview 
of the history of the problem in order to introduce the 
notation and to recall the underlying ideas. 

In 1883 Osborne Reynolds suggested the necessity of 
a statistical description of turbulence at high Reynolds 
number Re = VLjv, where V is a characteristic velocity, 
L is a characteristic length scale, and v is the kinematic 
molecular viscosity. Note that a typical value of Re for 
water flow in rivers and for air flow in the atmosphere is 
about 10 7 — 10 9 which is tremendously larger then the 
critical value Re cr ~ 10 2 at which a laminar flow loses its 
stability. 

The modern concept of hydrodynamic turbulence orig- 
inates from the Richardson cascade model (1922) Q] in 
which turbulent motions are produced because of the in- 
stability of the laminar flow around the streamlined body, 
the geometry of which determines the scale of the insta- 
bility L. The eddies with the characteristic scale L which 
appear due to the instability are also unstable in their 
turn and produce as a result eddies of smaller sizes, which 
are unstable again and produce eddies of even smaller 
sizes and so on and so forth. This process continues until 
the "current" Reynolds number dependent on the scale 
reaches the critical value Re cr . Eddies of smaller sizes 
are stable, and their energy dissipates into heat due to 
viscosity. In such a way at Re 3> Re cr fully developed 
turbulence arises in which there are turbulent motions at 
all scales from L down to some "viscous" scale r] <C L. 
We will refer to these motions at the scale r as r-eddies. 

Scaling behavior of velocity differences in fully de- 
veloped turbulence was predicted by Kolmogorov and 
Obukhov more than 50 years ago in their celebrated 
papers |^||. They assumed that the energy which is 
pumped on the largest scales ~ L is subsequently trans- 
ferred from larger eddies to smaller ones and that the 
statistics of the energy pumping is lost, except for the 
injection rate because it equals (in the stationary case) 
the value of the energy flux through each scale. In such a 
way the properties of turbulence in the inertial subrange 
of scales r, L r -q, are suggested to be universal 
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and to be independent both of details of the excitation 
of the turbulence and of boundary conditions. In partic- 
ular it was assumed that the statistics of fine scale turbu- 
lence (at r <C L) is homogeneous and isotropic. Clearly, 
the value of the energy flux is equal to the mean value 
e = (e(t, r)) of the energy dissipation rate e per unit mass 
occurring at small scales. 

Thus according to the Kolmogorov-Obukhov (1941) 
picture of fully developed, homogeneous, isotropic turbu- 
lence of an incompressible fluid (hereafter KO-41) there 
is only one relevant parameter e in the inertial subrange 
of scales which is the mean value of the energy dissipation 
rate. This means that in KO-41 a simultaneous correla- 
tion function of the velocity (and its derivatives) on a 
scale r in the inertial subrange is determined only by e 
and by r itself and is independent both of the outer scale 
L and of the viscous scale rj. This allows one to evaluate 
easily all correlation functions of turbulent fields in the 
inertial subrange with the help of dimensional reasoning 
(see for example textbooks 0k|). Some essential details 
are reproduced below, in Section |[ 

However the situation in the theory of turbulence is not 
as simple as it appear at first glance. In the 40's Lan- 
dau 0] indicated that the above KO-41 picture is not 
so obvious because of the intcrmittcncy inherent in tur- 
bulence. This means that there are relatively calm peri- 
ods which are interrupted irregularly by strong turbulent 
bursts either in space or in time. Experimental evidence 
of intermittency was given by Batchelor and Townsend 
in 1949 § and then by Kuo and Corrsin in 1971, 1972 
[0,1). Intermittency leads to strong fluctuations of the 
energy dissipation rate e, thus one may find amongst the 
parameters relevant for the description of turbulent mo- 
tions of scales r in the inertial subrange along with the 
mean value of the energy flux e, the dispersion of the en- 
ergy flux over the scale r. This dispersion may depend 
on the parameter L/r characterizing the number of suc- 
cessive crushings of r-eddies required to reduce the scale 
from L to r. 

A number of schemes were proposed to take intermit- 
tency into account. The first attempts belong to Kol- 
mogorov ||] and Obukhov [jlO| who proposed the so- 
called log-normal model of intermittency. After them 
Novikov and Steward [0|, Yaglo m Jl2| ], Mandelbrot 
[^3|,0, Frisch, Sulem and Nelkin [p^fand many oth- 
ers have proposed various modifications of the KO-41 
picture. These attempts to construct phenomenologi- 
cal models are very important as a method to represent 
experimental data and interesting from the theoretical 
viewpoint. However all of them suffer from a lack of con- 
nection to the Navier-Stokes equation which describes 
the real fluid dynamics. 

Various closure procedures have been also constructed, 
see for example Obukhov 1941 @, Heisenberg 1948 [p7| , 
Kraichnan 1965 p8[ , etc.; for a review see, i.e., |l{|. Af- 
ter the exclusion of non-physical approximations one can 
claim that all these theories give KO-41 scaling of ve- 
locity differences. However all such closures are uncon- 



trolled approximations since there is no small parameter 
in the theory of turbulence. Therefore such approxima- 
tions may be considered as a tool to evaluate some dimen- 
sionless constants, but not as a proof of KO-41. Indeed, 
one can say that intermittency corrections possibly stem 
from some spatial structures of turbulent flow, like fila- 
ments, etc. The existence of such structures is reflected 
in the high moments of velocity differences and therefore 
cannot be described within the framework of any closure 
procedures. Attempts to reach this phenomenon may be 
done in a systematic theory of turbulence without any 
truncations. 

The first attempt to formulate a systematic statistical 
description of turbulence directly from the Navier-Stokes 
equation was made by Wyld |20(| . Another derivation 
of the very same diagrammatic approach was suggested 
by Martin, Siggia, Rose J21J in the framework of a func- 
tional integration approach. Unfortunately these tech- 
niques lead to the representation of the correlation func- 
tions in the form of an infinite series and any truncation 
of these series breaks the Galilean symmetry of the prob- 
lem. It leads to enormous technical difficulties. Formally 
these arise from the infrared behavior of integrals which 
become more and more divergent with increasing order 
in the perturbation series. The physical reason for this 
is the sweeping effect of eddies of given scale in the in- 
ertial interval by the velocity field of all eddies of larger 
scales. In order to eliminate this sweeping effect from 
the theory Kraichnan formulated p2[ a perturbation ex- 
pansion in the terms of the Lagrangian velocity. Result- 
ing Kraichnan's "Lagrangian-history" renormalized ex- 
pansion is a systematic perturbation procedure which 
possesses Galilean symmetry order by order. In this 
approach the "Lagrangian-history direct-interaction" ap- 
proximation is just the first step p8| . There is the price 
to pay for the elimination of the sweeping in this the- 
ory: The "Lagrangian-history" perturbation expansion 
has NO Feynman-type diagrammatic representation of 
n-th order terms in the expansion. Therefore it is very 
difficult to analyze high-order terms in this approach and 
we do not know any example of its use in a theory of 
hydrodynamic turbulence besides direct interaction ap- 
proximations. 

A different way to overcome the problem of infrared 
divergences due to the sweeping is by the use of the 
Belinicher-L'vov resummation p3[ | of the Wyld diagram- 
matic expansion. This resummation corresponds to a 
change of variables in the effective action of the Martin- 
Siggia-Rose approach |^l[ from the Eulerian velocity to 
the quasi-Lagrangian (qL) velocity. The qL description 
of turbulence was suggested by L'vov in 1980 (sec p3fl ). 
The relation between the qL velocity u and the Eulerian 
velocity v is 

v(t,r) = u(t,r- e (t)) (0.1) 

where the time evolution of the vector g(t) is determined 
by the equation 
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dg/dt — u(t, r ) 



(0.2) 



where ro is a marked point which may be arbitrarily cho- 
sen. Let us stress again that there are no approximations 
in this step; the above relation is simply a change of vari- 
ables in the Martin-Siggia-Rose effective action. A sys- 
tematic analysis of resulting resummed series was done 
in paper p3| ]. In contrast to Kraichnan's Lagrangian- 
history approach the resulting series do have a Feynman- 
type diagrammatic representation. In contrast to the ini- 
tial Wyld diagrammatic technique |2(J (and the Martin- 
Siggia-Rosc one |2l|]) resummed diagrammatic series do 
possess Galilean invariance order by order before and af- 
ter Dyson's line renormalization. There is a price to pay 
for these advantages: the transformation ( |0.l| , |0.2| ) breaks 
the space homogeneity of the problem. Therefore the 
Green's and correlation functions in the qL representa- 
tion depend on the space coordinates r and r' separately 
in contrast to standard approaches in which these func- 
tions depend only on their differences. However it was 
possible to overcome the corresponding technicalities |2|] 
and to prove that there are no infrared nor ultraviolet di- 
vergences in any term in any order of perturbation series. 
As a result we have a powerful technique which allows one 
to reach, as we believe, crucial progress in understanding 
the scaling of developed hydrodynamic turbulence. 

The important step on that way was made in |]23|| , (see 
also |24|]) where it was shown that the structure func- 
tions of velocity differences with K041 scaling are an 
order by order solution of the corresponding diagram- 
matic equations. In our 1993 paper |2j| we show that 
the KO-41 scaling of velocity differences is the unique 
solution (in some region of scaling exponents) under the 
very plausible assumption that the time-dependent corre- 
lation functions of qL velocities are scale invariant. This 
is a non-perturbative result which follows from our "fre- 
quency sum rule for the dressed vertex" ^5|. The later 
is an exact consequence from the causality principle for 
the three-particle Green's function. 

If we do believe that the Belinicher-L'vov resumma- 
tion [ p3|]2~i| ] gives an adequate formalism to describe de- 
veloped hydrodynamic turbulence we should next dis- 
cover in these terms the mechanism for anomalous (non- 
Kolmogorov) scaling of the energy dissipation field. We 
have introduced such a mechanism in brief p6| . In 
the direct qL diagrammatic expansion for the correla- 
tion function of the energy dissipation K £e (R) (1.6) we 
found an infinite subset of logarithmically diverging dia- 
grams. Their resummation leads to a renormalization of 
the normal Kolmogorov-41 scaling behavior of K ee (R). 
Let us stress that this is again a non-perturbative result 
obtained with a diagrammatic technique. There is no 
closure procedure which leads to anomalous scaling; and 
only resummation of an infinite series of relevant terms 
of the diagrammatic perturbation expansion may pro- 
duce this effect pq |. In "physical language" the mech- 
anism responsible for the anomalous scaling was called 
the telescopic multi-step eddy interaction. The essence 



of this mechanism is the existence of very large number 
(R/n) Aj ^ 1 of channels of interaction of large eddies of 
scale R in the inertial interval with eddies of viscous scale 
rj via a set of eddies of all intermediate scales between R 
and rj. Note that the same mechanism works also for any 
local field constructed from the velocity gradients, e.g., 
w 2 where u = V x v is the vorticity vector. Our find- 
ings mean that KO-41 scaling describes only a part of a 
full set of correlation functions characterizing developed 
turbulence. 

The next challenge for the resummed description of 
turbulence is to understand the observable deviations of 
the scaling of structure functions of velocity differences 
from the KO-41 prediction. Although a consistent the- 
ory of this phenomenon belongs to the future a possible 
way to describe these deviations as intermediate asymp- 
totic behavior structure functions at very large but finite 
Reynolds numbers is suggested by L'vov and Procaccia 
in their "subcritical scenario" [^7) . If we accept this sce- 
nario it would mean that all available experimental obser- 
vations of the scaling behavior of developed turbulence 
may be understood in the framework of the resummed 
diagrammatic approach p3j. 

In the present paper a consistent theory of the anoma- 
lous scaling of the energy dissipation field (and some 
other hydrodynamic fields constructed from velocity 
derivatives) is developed in the limit Re— > oo. The the- 
ory begin with the Martin-Siggia-Rose approach (2l|]) to 
the Navier Stokes equation and based on the Belinicher- 
L'vov resummation of corresponding diagrammatic se- 
ries. It gives an analytical description of the telescopic 
multi-step eddy interaction suggested in our Letter p6| . 
We show that the correlation functions of the hydrody- 
namic fields constructed from velocity gradients possess 
a complicated behavior characterized by an infinite set 
of scaling dimensions which do not reduce to the Kol- 
mogorov ones. Unfortunately these exponents cannot 
be calculated explicitly since they describe the scaling 
behavior of a complicated integral equation the kernel 
of which is not found analytically. Nevertheless a set 
of relations the between scaling behaviors of correlation 
functions of different hydrodynamic fields can be estab- 
lished. We found a number of selection rules related to 
the symmetries, including also the hypothetical confor- 
mal symmetry of developed turbulence, and formulated 
some restrictions related to the incompressibility condi- 
tion of hydrodynamic motion. Our scheme leads also to 
an anomalous asymptotic behavior of correlation func- 
tions of velocity differences which are described by the 
same scaling exponents. 

A tool which is useful for us is the set of so-called fu- 
sion rules for fluctuating fields introduced by Polyakov 
p8[ in the context of second order phase transition the- 
ory. These rules enable one to determine the asymptotic 
behavior of a correlation function in the case where a 
number of points (in the real r-space) are separated from 
other points by a large distance. Kadanoff [^9| and Wil- 
son fpol formulated these rules in the form of the so-called 
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operator algebra. This algebra was extensively used in the 
conformal theory of 2d second order phase transitions 
fjjj . Since the algebra is not related to any particular 
properties of the system one may expect that it can be 
formulated for any system described by a set of corre- 
lation functions with a scaling behavior. In this work 
we are going to propose an operator algebra for 3d tur- 
bulence, which is characterized by a set of correlation 
functions of fluctuating turbulent fields. 

The structure of the paper is as follows. Section | 
contains a brief overview of the basic physical concepts 
of hydrodynamic turbulence associated with the idea of 
scaling. In Section [n] we develop the qualitative picture 
of the "telescopic multi-step mechanism" of eddy inter- 
action suggested in our 1994 Letter j26). This allows one 
to recognize physical reasons for the anomalous scaling 
behavior of the energy dissipation field and other hydro- 
dynamic fields. In Section III we develop the analyti- 



cal theory beginning with the Navier-Stokes equation in 
the quasi-Lagrangian representation and based on a line- 
renormalized diagrammatic approach. We extract infi- 
nite sub-sequences of diagrams with ultraviolet logarith- 
mic divergences which lead after summation to anoma- 
lous scaling behavior of the correlation functions of the 
hydrodynamic fields with different scaling exponents. In 
Section fV we formulate rules, including the fusion rules, 
which allow us to predict scaling behavior of correlation 
functions of hydrodynamic fields constructed from the 
velocity gradients and of the velocity differences. The re- 
sults obtained and a discussion concerning the region of 
applicability of our approach are presented in the Con- 
clusion. 



I. SCALING IN HYDRODYNAMIC 
TURBULENCE 

In this section we recall the main ideas associated with 
the scaling behavior of different correlation functions of 
the hydrodynamic fields such as velocity differences and 
fields constructed from the velocity gradients. We in- 
troduce the notion of normal scaling behavior determin- 
ing by the Kolmogorov's dimension estimates and also 
anomalous exponents which describe deviations from this 
normal behavior. 



A. Kolmogorov 1941 Scaling for Turbulent Velocity 

Let us begin by accepting The Kolmogorov's 1941 
picture of turbulence (KO-41) and consider within this 
framework the statistical properties of the turbulent ve- 
locity field v a (t , r) . The principal quantity characterizing 
the turbulence in KO-41 is the average value e of the en- 
ergy dissipation rate 



where s 2 is the second power s 2 = s a ps[j a of the strain 
tensor. s a p — (y a vp + V^i> Q )/2. Note that e is the 
energy dissipation rate per unit mass; to find the energy 
dissipation rate per unit volume e should be multiplied 
by the mass density p. 

We should be careful in treating correlation func- 
tions containing the velocity field v(t, r) itself. First, 
this quantity is not Galilean invariant and depends on 
the choice of the reference system. Second, the main 
contribution to correlation functions of velocities (like 
(v(i, r) • v(t, r'))) comes from the largest eddies. This 
means that expressions for such correlation functions 
contain large homogeneous contributions which cannot 
be obtained in terms of the fields defined in the iner- 
tial interval. As a result such correlation functions are 
not universal: they may depend on boundary conditions, 
statistics of pumping etc. In the analytical approach the 
dominating contribution to correlation functions associ- 
ated with the energy-containing subrange is reflected in 
infrared divergences of corresponding integrals. It is well 
known that to avoid this difficulty we should consider 
such objects as velocity differences which are Galilean 
invariant; their correlation functions do not contain in- 
frared terms related to the largest eddies. 

Let us introduce the traditional designation for the av- 
erage values of the powers of the velocity difference: 



D n (R) = (\v(t,r)-v(t,r + R) |") . 



(1.2) 



The quantities D n are usually called the structure func- 
tions of the velocity differences. We will be interested 
in the behavior of these structure functions in the iner- 
tial subrange L ^> R ^> r\. In Kolmogorov's picture one 
can evaluate D n (R) as e a R b with the exponents a and 
b chosen in such a way as to have the same dimension- 
ality (remind that s is the energy dissipation rate per 
unit mass). This condition gives a = b = n/3, which 
leads to the famous KO-41 result D n (R) = C n {eR) n / 3 
with some dimensionless constants C n . Physically these 
dimensional estimates mean that the principal contribu- 
tion to the functions D n (R) is associated with i?-eddies 
having characteristic velocity V R of the order of 



Vr, 



(sR) 



1/3 



(1.3) 



The characteristic turnover time of i?-eddies in KO-41 
can be estimated as R/V R . Since in KO-41 one can con- 
struct only one quantity with the dimensionality of time, 
the lifetime of i?-eddies t r (at R in the inertial subrange) 
should be of the order of their turnover time. 

In treating correlation functions of the velocity gra- 
dients one has no problems with large-scale terms since 
these correlation functions are not sensitive to motions 
on the pumping scale L. KO-41 dimensional reasoning 
allows one to find the scaling behavior of the two-point 
correlation function of the strain tensor s a p and the vor- 



e(t,r) = 2us 2 (t,r) , 



(1.1) 



ticity vector u> = V x v. For this we can use (1.3) which 
gives for eddies of the scale R that s,u> ~ e i/a i? _2//3 . 
Moreover, one can establish the tensor structure of the 
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pair simultaneous correlation functions of s,w starting 
from 



D a0 (R) = (MR) - u a (0))MR) - v p (p))) 
= ^C 2 sR 2/3 (6 a(j - R a Rp/4R 2 ) . 



(1.4) 



The tensorial structure here is determined by the incom- 
pressibility condition V • v = 0. C'2 is a dimensionless 
constant. The pair correlation functions of the velocity 
differences can now be found as 

(V M »„(R)V^(0)) = -V^A^R) . (1.5) 



B. Normal and Anomalous Scaling in Turbulence 

In the previous Section the "naive" KO-41 scaling ex- 
ponent for some correlation functions have been intro- 
duced. We refer to these exponents as normal. As we 
noted strong fluctuations could change these "naive" ex- 
ponents. We refer to the deviations from KO-41 values 
as anomalous exponents. Accordingly we introduce in 
this Section the formal designations for the anomalous 
exponents characterizing different correlation functions. 

Intermittency is associated with strong fluctuations of 
the energy dissipation rate e. To characterize these fluc- 
tuations we should treat correlation functions of s. In- 
troduce the designation for the pair correlation function 



K ee (R) = ((e(t,r)e(t,r + R))} 



(1.6) 



Here and henceforth the double angular brackets mean 
the irreducible part of the correlation function of two 
fields: 

(ty(nMr 2 )» = (V(riMra)) - (V>><¥>> • 

It is natural to assume [0 that the correlation function 
K ee (R) is scale invariant in the inertial subrange: 



K es (R) oc R-* 1 



(1.7) 



with some phenomenological exponent /i. This quantity 
is usually called the intermittency exponent. 

There are two known ways to evaluate the exponent 
/i within dimensional reasoning. The simplest one is as 
follows. Clearly the dimensionality of K £e (R) coincides 
with the dimensionality of e 2 . Therefore an attempt to 
express K ££ via powers of e and R gives p = 



K se (R) ~ e 2 R° 



(1.8) 



In contradiction with the folklore (see, for example ||) 
thi s est imate is not related to KO-41. Indeed, according 
to (1.1) the energy dissipation field s(t,r) is proportional 
to the viscosity v which is not a value characterizing the 
inertial subrange and it is not reasonable to rewrite the 



parameter v characterizing the medium in terms of s and 
r. To be consistent with KO-41 we should treat the corre- 
lation functions independently of v, in terms of s 2 . Then 
according to KO-41 dimensional reasoning 



«s 2 (r) S 2 (r + R)}} ~g*/ 3 ir 8 / 3 . 



(1.9) 



ws"s~ this expression means that fi — 8/3. 
The explicit expression for K ee (R) with /j, = 8/3 was 
found by Golithin |52| under the assumption of Gaus- 
sianity of the velocity fluctuations. As we show one can 
obtain it within the "spirit" of KO-41. We see now that 
the normal KO-41 exponent leads to an extremely fast 
decay of se correlations. 

During last three decades a lot of effort has been ex- 
pended in attempting to measure the exponent n (see 
e.g. [^3|-^ll). There are many experimental difficulties in 
this problem, therefore different authors give for /1 val- 
ues ranging from 0.1 to 0.5. Thus it is clear that the 
normal KO-41 exponent \i = 8/3 does not correspond to 
experiment. This is not amazing at all, because from a 
consistent theor etic al viewpoint there is no re ason to ac- 
cept both Eq. ( [Tg| ) giving fi — 8/3 and Eq. (1.8) giving 



fx = 0. Note that the multifractal pictures of turbulence 
with some assumption about the structure of the flow 
|^2| , ^3| give for the exponent p an acceptable value nea r 
1/3. We conclude that the correlation functions (1.6) 
undoubtedly possess an anomalous behavior. 

The physical rea son for deviations of the correlator 
K es (R) from (l.£) is clear, see for example ^JDsf . 
Namely, the s 2 -, and the w 2 -fields as well as s 2n -ficlds 
and many others, are so-called viscous subrange fields in 
the sense that mean values of these fields are mainly de- 
termined by eddies of characteristic scale 77. It means 
that the 77-eddies may play an essential role in forming 
the behavior of the correlation functions of these fields. 
Therefore the correlation functions of these fields have 
to be not only a function of the inertial subrange val- 
ues as in the normal KO-41 estimate (l.E), but also a 
function of the viscous scale rj. As a result powers of 
the dimensionless parameter R/r/ can enter as a factor in 
all correlation functions of s 2n and s 2ra -fields. There is 
no dimensional reasoning which can determine to what 
power this parameter appears. So, the only result which 
follows from a consistent dimensional reasoning is 



((s 2n (K)s 2m (0))) = C n 



\R 2 



2(m+n)/3 



(1.10) 



with some unkn own expon ent s A nm . Note that both 
particular cases (U) and (|J) are described by ( |TTTo| ) 
f or so me values of An. We will consider the relations 
( 1.10 ) as a definition of the anomalous scaling expo- 
nents A nm which determine the scaling behavior of the 
prefactor of the naive KO-41 expression for the correla- 
tion function. Exactly the same problem with unknown 
dimensionless factors arises for many-point correlation 
functions of other viscous subrange fields. An impor- 
tant question concerns also the dimensionless factors rj/R 
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and R/L in the structure functions of velocity differences 
(1.2). The appearance of such factors is not forbidden 
from the general point of view, if to i nclu de these factors 
besides KO-41 estimates we find for (1.2): 



D n {R) oc R 1/3 



(1.11) 



where the scaling exponents Cn differ from the KO-41 
prediction n/3. In his lognormal model of intermittency 
Pj Kolmogorov argued that 



Cn = n/3 — fj,n(n — 3)/ 18 



(1.12) 



In the Novikov-Steward |llj model and in the /3-model 
by Frisch, Sulem, and Nelkin 



Cn = n/3 - fi(n - 3)/3 



(1.13) 



Here \i is the same scaling exponent as in K e£ (R) oc R~^. 
Note that in these models the deviation of Cn from n/3 
is related to the decay of ee correlations. In fact, there 
is no solid basis for this statement ]46| , f47| ] . 

Actually we do not see serious theoretical reasons for 
the deviations of Cn from their KO-41 values n/3 in the 
limit Re— > oo. Indeed, on one hand there is the theo- 
rem proved by Belinicher and L'vov |23|,^j who demon- 
strated that the KO-41 structure functions of velocity 
differences order by order satisfy the corresponding dia- 
grammatic equations. Moreover, in 1993 p5fl we showed 
that KO-41 is the unique solution (in some region of 
scaling exponents) under the assumption that the time- 
dependent correlation functions of qL velocities are scale 
invariant. On the other hand we do not know any mecha- 
nism for the renormalization of the scaling in the absence 
of ultra-violet and infrared divergences of diagrams in 
nonlinear problems with strong interaction. Therefore in 
the present paper we have accepted the KO-41 scaling 
of velocity differences (Cn = n/3) and find some remark- 
able consequences of this fact concerning the anomalous 
exponent of different correlation func tions of the veloc- 
ity gradients as introduced in ( 1.1C ). Certainly there 
are deviations from Cn = n/3 which are observed experi- 
mentally (see, e.g., ]48|,[49| and references therein) and in 
numerically pC[ |. However we believe that these devia- 
tions are related to the finite value of Re in experiments 
which leads to restricted values of the inertial subrange 
L/rj < 10 4 for the largest available Reynolds numbers 
Re < 10 s — 10 9 (see also Conclusion). A possible "subcrit- 
ical" mechanism for such intermediate non-Kolmogorov 
behavior was suggested recently by L'vov and Procaccia 



II. TELESCOPIC MULTI-STEP EDDY 
INTERACTION. 

Our aim in this Section is to evaluate the correlation 
functions of s 2 and J 1 for a separation distance R within 



the inertial subrange L ;§> R ^> n. We show that the lead- 
ing contributions to the correlation functions K S 2, S 2(R), 
K S 2,u2{R), and K U 2,ui{R) have anomalous scaling be- 
havior characterized by the same scaling exponent Ai. 
We argue that at the same time the correlation function 
of a traceless tensor (say uJaiop — ^5 a 0U! 2 ) has to have 
another independent anomalous scaling exponent. The 
telescopic multi-step eddy interaction mechanism for the 
anomalous exponent of the energy dissipation field was 
suggested in our Letter (2(| through the diagrammatic 
approach. In order to elucidate the physical basis of the 
involved diagrammatic calculation we first describe our 
findings in this Section using the popular handwaving 
language of cascades, eddies and their interactions before 
turning to the cumbersome technicalities of the analyti- 
cal theory of turbulence. 

In our approach K a 2 >S 2 and other correlation functions 
are represented as a sum over the contributions from tur- 
bulent fluctuations in the inertial range. The fluctuations 
s 2 (r) and uj 2 (r) are due to fluctuations of the velocity gra- 
dient at point r which in turn may be regarded as the 
result of a superposition of shear rates stemming from 
velocity differences on various length scales r, which ac- 
cording to KO-41 scale like r~ 2 / 3 . Therefore the main 
contribution to the fields of interest themselves is ex- 
pected to come from the smallest scales near rj while the 
main contribution to the correlation over the distance R 
is naively expected from the smallest scales which bridge 
the gap between r and r + R, namely from the scales R. 
This means that scales larger and smaller than R can be 
neglected in our considerations. We call the above rea- 
soning naive because it tacitly assumes that an indirect 
action through scales < R is much weaker. 

In this Section first we consider the simplest contribu- 
tion to the correlation function of our fields coming from 
eddies of one scale R. Than we show that this is not 
the whole story and that the combined effect of indirect 
interactions of scale R with scale n is also important. We 
will examine the contribution from two and three groups 
of eddies of different scales. Finally we consider the total 
contribution of eddies of all scales from n to R and show 
that the indirect interactions of scale R with scale r\ via 
the set of all intermediate scales gives rise to a correction 
of the scaling exponents of our fields from Kolmogorov's 
values downward. 



A. Direct Contribution of One-Scale Eddies 

In order to formally determine the notion of x-eddies 
let us partition fc-space into shells separated by wavenum- 
bers k n = 2it/x n . By n-eddies (or eddies of scale x n ) we 
mean the turbulent velocity field v n (t, r) which has only 
part of the Fourier harmonics of the "real" turbulent ve- 
locity field v(£, k) with k between k n -i and k n . 

Consider now the contribution to K S 2 tS 2(R) due to ed- 
dies of some characteristic scale x„ . It is clear that for a 



G 



separation distance R which is about the eddy scale x n , 
the correlation functi on K S 2, S 2(R) has to be of order of 
the KO-41 estimate ( |l.9| ). Indeed, in this case R and e 
are the only parameters in the problem considered. 

In the case of x n 3> R (but still x n <C L) the correlator 
Ks2,s2(R) is independent of R and is approximated by 



K 



s2,s2 



(X n ) 



1/3 



Thus this estimate of 



Ks2,s2{R) is smaller than ( |l.9| ) 

For a;„ <C R the correlator ifg 2.s2( -R) has to be ex- 
ponentially small with respect to (1.9). Therefore for a 
one-scale contribution the largest estimate of K ee {R) is 
give n for x n — R and corresponds to the KO-41 value 
& 



B. Contribution of One-Step Eddy Interactions 

Consider now the contributions of two group of ed- 
dies (or scales x n and x m ) homogeneously distributed in 
space. First of all there are the contributions of each 
group of eddies separately. Second there is a cross- 
contribution arising from the interaction of these two 
groups of eddies. By discussions similar to those pre- 
sented in the previous Subsection one may see that the 
largest cross-contribution to K s2 , S 2(R) in this case ap- 
pears at x m — R, and an arbitrary value of x n (or vice 
versa). Therefore we choose from the beginning m = N 
with an arbitrary value of n. 

Let us consider correlations between fluctuations on 
different scales with the help of a conditional proba- 
bility density. The larger eddies modify the statistics 
of the smaller ones locally in r-space by their rate- 
of-strain field. Thus the probability density to find 
on scales < n near r the velocity gradient Vv is 
P„(Vv(r + r')|V(v„ + i(r)+v„ +2 (r) + ...)), where the 
condition is set by all the larger scales and jr'| <C x n +\. 
We assume that fluctuations around the local mean value 
at points farther apart than |r'| <C £ n +i are not corre- 
lated. Further w iU denote such conditional av- 
erages. Obviously, |_L---J' m _|„ = L---J( holds with I = 
max(m, n). KO-41 scaling implies that the velocity gra- 
dients of larger scales are comparatively small. Hence, 
we will make use of the expansion of P n for small large- 
scale gradients, P„(Vv|Vv >n ) = Pi 0) (Vv) + P^Vv) • 
Vv>„ + Pi 2) (Vv) • Vv > „Vv > „. 

We illustrate in this subsection the consequences of 
such correlations considering the hypothetical case where 
only fluctuations in shell N and in shell n, 1 <C n <C N 
are excited. Then 



K 



s2,s2 



(\r 1 -r 2 \) = ((s 2 (t,v 1 )s 2 (t,v 2 ))) 



s 2 (r 1 ,t)s 2 r 2 ,t) 



L^(*,ri)J B LP(t,r a )J B J 



(2.1) 



where we denote s 2 (i,r) — (s 2 ) as s 2 (i,r). The last step 
in (|2.lD is justified because the s 2 field is mostly sensitive 



to the fast, small scale motions; each x-ensemble aver- 
age can be done in the presence of some realization of 
the P-eddies. and only when we compute the correlation 
function (2.1) we need to average in the P-ensemble. In 
other words one may split the averaging in the above 
equation because the correlation length under the condi- 
tional n-average is much shorter than R. 

Consider now the average over small n-eddies 



(2.2) 



Using the expansion of the conditional probability we 
obtain for this object an expression of the form 

[?(t, r)J n ~ si (t, r) + A a ^ n V a v 0)N (t, r) 

+B a ^s,n V a wft s (t, r) V 7 W5,„ (*, r) (2.3) 
-|-higher order terms . 

The first term on the RHS of this equation is the di- 
rect contribution of the TV-shell to Le(r)Jn- The largest 
contribution to [ s2 ( r )Jn comes from the n-shell itself, 
and is |_Sn( r )J™- However this contribution is indepen- 
dent of time and space coordinates and is canceled in 
the subtraction of e. The last two terms derive immedi- 
ately from the expansion of P n . The term of first order in 
Vv is not important because all contributions originating 
from this term will vanish finally under the average over 
fluctuations on intermediate scales. The expansion coef- 
ficients B are the corresponding derivatives of the func- 
tion [s (t, r)\ x taken at zero value of v N (t, r). Therefore 
they reflect the "intrinsic" properties of n-scale turbu- 
lence, without interacti on w ith i?-eddics (that is the rea- 
son why we display in ( |2.3| ) the dependence of A and B 
on the scale x n ). Note that the turbulence of x-eddies 
is isotropic and because of this the matrices A and B 
have to have some particular form which allows one to 
represent the above expansion in a more elegant form: 



b 2 (i,r)J, 



f Bn,„5(*> r ) + B\2,nU^ N (t, r) + ... (2.4) 



where s 2 N (t, r) and uu 2 N (t, r) are the contributions to the 
square of the strain tensor and the vorticity coming from 
the TV-shell and " " denotes the irreducible part in 



accordance with (2.2). It is very easy to understand why 
one ha s only two terms (oc s 2 and cx to 2 ) in the RHS 
of (2.4): the LHS of this equation is scalar (under all 
rotations) therefore the RHS also has to have a scalar 
form, and there are two scalars (s 2 and to 2 ) which one 
can build from (V a! W/3)(V 7 u ( 5). Clearly, the exp ansi on 
for another scalar \uj 2 (t,r)\ x has to be similar to (2.4): 



L„ 2 (t,r)J n =„ 2 (M) 

+B 2hn s%(t,r) + B 22>n i? N (t,r) + ... . (2.5) 

One may diagonalize the matrix Bij and find linear com- 
binations of fields s 2 and lo 2 
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*i(t,r) = U n s 2 + U 12 u;^ 
* 2 (t,r) = U 21 s 2 + U^uj 2 



(2.6) 



for which the expansions ([Q.2.5) (up to second order 
terms in Vv) take the simplest form: 



L*i(*> r )J„ 

* J - iV (t J r)[l + S i , n ] 



(2.7) 



J 



1,2 . 



The coefficients Bj n can be estimated applying physi- 
cal reasoning (or, equvalently, the Navicr Stokes equa- 
tions). The dimensionless parameter describing the rel- 
ative change of velocity v„(r) with varying Vv w is 
T n Vv N (r), where r n is the life time of n-eddies. There- 
fore, Bj. n ~ r 2 |_|Vu n | 2 Jn ~ %„■ Here we have used essen- 
tially the KO-41 result that the life time r n is approxi- 
mately the turnover time of the n-eddies. This allows one 
to conclude that Bj jTl is independent of the scale x n . This 
statement is of crucial importance for the understanding 
of the origin of the anomalous scaling: all scales in the 
interval from R to r\ contribute equally to the correla- 
tion functions under consideration. Therefore we have to 
take into account not only contributions from two group 
of eddies, as we did up to now, but contributions of all 
the group of eddies of scales from the above interval. 



jv 



(2.10) 



where N = log A (i?/?y), Aj = ln(l + 5,-)/ln(A) and 
j = 1, 2. Following the terminology of second order 
phase transitions one can call the exponents Ai and A 2 
anomalous dimensions of the fields an d ^ 2 . The 
physical meaning of Aj is clear from ( ^.10| ): the factor 

(i?/f?) Aj is the total number of effective channels of multi- 
step eddy interactions of the ^f(R\t, r) field (having scale 
R) with the smallest 77-eddies. 

Return ing t o the original fields s 2 and oj 2 one has in- 
stead of fl2.10j ): 



'R\ A i 



~ A n V 1>N {t,r)(-) + A 12 # 2 ,„2(i,r)(-J , 

l...^(t,r)\ y \ x ...\ (2.11) 

~ A a i*i )W (t, r) (-J + A 22 V 2<N (i, r) (-) 
with some dimensionless coefficients An . 



C. Contributions of Two- and Many-Step Eddy 
Interactions 



D. Additional Anomalous Fields and Correlation 
Functions 



Now let us consider the contribution of three groups 
of eddies (with scales R, x n , and x m ). This will directly 
lead to anomalous scaling. For this contribution one has: 



L*i(*,r)J„ 



(2.8) 



The first three terms on the RHS of this equation describe 
the direct contribution of the i?-shell and the contribu- 
tions from the direct interaction of the i?-shell with the n- 
and m-shell, respectively. The last term (oc B^ n Bj_ m ) is 
due to the indirect effect of the largest scale, the -R-shell, 
on the smallest scale, the n-shell say, via the intermedi- 
ate TO-shell. To obtain this term one has to repeat twice 
the above expansion. The RHS of fl2.8| ) is proportional 
to (1 + Bj, n )(l + Bj jn ). Thus it is plausible that one 
obtains in the case of N shells the result 



JV 



L*j(*,r)L 



J,N 



(*>*•) n 



1 + B 



(2.9) 



Note that the independence of the Bj^ n of n as pointed 
out above will only hold when the width Ak n of the shells 
scales in the same way as the k n themselves. We choose 
k n +i/k n — A so that neighboring shells may be consid- 
ered as almost statistically independent. Such A > 1 
does exist because of the locality o f ene rgy transfer via 
the scales p3i . Then one can write (|2.9|) in the form 



Consider now correlation functions of the pseudovector 
field f*= SapoJp- Cle arly, one can expand [f*(t,r)\ x in 
a similar way to (2.4). Previously we expanded a scalar 
in terms of scalars. The pseudovector field has to be ex- 
panded in the terms of pseudovector fields. Generally we 
have to expand a field with given transformation proper- 
ties (under all rotations and inversion) in terms of fields 
with the same transformation properties. This require- 
ment is a consequence of the isotropy of turbulence and of 
the fact that the expansion coefficients reflect the prop- 
erties of the fine scale turbulence only. So, 



Lf*(i,r)j*~f*(i,r)[l + &* + 



(2.12) 



Here b* is some new dimensionless expansion coefficient. 
Repetition of all of above considerations allows one to 
conclude that 



[...[r(t,r)\ y \ x ...\^r N (t,r)(^ A 



(2.13) 



with new anomalous scaling exponents A* oc ln(l + b*). 
In the same way one can find three new anomalous scaling 
exponents, associated with three traceless tensor fields 
of the second order (in Vv), one new exponent for a 
pseudotensor of the third order and so on. For more 
detail, see Section IVB. 

Now one can proceed to establish the scaling behavior 
of the two-point correlation functions of the s 2 -, w 2 -, /*-, 
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etc. fields. For example, Eqs. (2.11) together with (2.1) 
allow one to see that 



K 



«2.«2 



(R) 



g 4 / 3 

R 8 / 3 



i?\2Ai 



(2.14) 



i?\ 2A 2 



F 4/3 



1 

A A fR\ 2Al 



(2.15) 



R\A X +A 2 /i?\2A 2 

+ A12A22 — 
T] / V 77 / 



and a similar equation for K U 2,u2- The leading term in 
each of th ese expressions scales like i? 2Al ~ 8 / 3 . Expres- 
sion (2.11) also allows one to estimate in the same manner 
cross-correlation functions of s 2 and to 2 with other hydro- 
dynamic fields, to estimate the fjT^ -correlation function 



f *,f * 



«f*(ri)f*(r 2 )» 




+ ... , (2.16) 



etc. We will discuss the problem of the evaluation of 
the correlation function of various hydrodynamic fields 
in more detail in Section |y|. 



E. Summary of the "Physical" Approach 

One should to realize that the consideration of anoma- 
lous scaling in this Section is just an attempt to explain in 
terms of cascades, eddies and their interaction the com- 
plicated multi-step processes in a system with strong in- 
teractions without a small parameter. It has not been 
a simple task because our "physical" intuition is usually 
restricted to the first or second order of the perturba- 
tion approach while the origin of strong renormalization 
of scaling behavior is related to taking into account an 
infinite series of terms in perturbation theory. In result 
our "physical" approach of this Section cannot be consid- 
ered as consistent. For example in the beginning of our 
"explanation" we discussed some group of eddies with 
well separated scales R ^> x ^> y 3> 77. Then it was 
allowed to average our fields in a few consequent steps, 
over fy-ensemble, then over y-, x-, and finally over i?- 
ensemble and to make use of expansions of the fields in 
powers of the gradients of the velocity fields, taking into 
account just the first non-trivial term. Then the main 
contribution to the renormalization of the scaling came 
from the multi-step processes in which our " group of ed- 
dies" are packaged densely in the inertial range of scales. 
Therefore the applicability parameter of this "physical 
way of thinking" which should to be small is really of 
the order of unity. This does not change qualitatively 
the physical picture of the anomalous behavior but leads 



to some additional contribution to the correlation func- 
tions which may be important. One direct way to find all 
of these contributions and to develop consistent descrip- 
tion of anomalous scaling in hydrodynamic turbulence is 
to make use of the diagrammatic perturbation approach. 
That is the topic of the next Section. 



III. DIAGRAMMATIC ANALYSIS 

The theory of turbulence is a theory of strongly fluctu- 
ating hydrodynamic motions. Systems with strong fluc- 
tuations are examined in quantum field theory (where 
quantum fluctuations are relevant) and also in condensed 
matter physics, e.g., in treating second order phase tran- 
sitions (where classical thermal fluctuations are relevant). 
Note that in both quantum field theory and the theory of 
second order phase transitions a scaling behavior of cor- 
relation functions (like in turbulence) is observed. It is 
known from these theories that adequate tools of theoret- 
ical investigation of strong fluctuating systems are based 
upon functional integration (path integration) methods, 
on different versions of the diagrammatic technique and 
on related methods like renormalization-group analysis. 
Therefore a consistent theory of turbulence should also 
be constructed in these terms. In this section we demon- 
strate how the diagrammatic analysis leads to the ap- 
pearance of anomalous dimensions in the theory of tur- 
bulence. 



A. Scaling and Divergences 

The diagrammatic technique for the problem of tur- 
bulence was first developed by Wyld p0[ |, who started 
from the Navier-Stokes equation with a pumping force 
on the right-hand side. The Wyld technique enables one 
to represent any correlation function characterizing the 
turbulent flow as a series over the nonlinear interaction. 
Unfortunately infrared divergences appear in the tech- 
nique. Physically they are related to the sweeping of 
small eddies by the velocity of the largest eddies. How- 
ever the sweeping does not change the energy of small ed- 
dies and consequently does not influence the energy cas- 
cade. This means that the form of the simultaneous cor- 
relation functions of the velocity is not sensitive to sweep- 
ing. In order to see this fact in each order of diagram- 
matic expansion we make use of the Belinicher-L'vov re- 
summation p3| which corresponds to the strongly non- 
linear change of variables: from the Eulcrian velocity to 
the quasi-Lagrangian ones ([Xl.0.2). It was demonstrated 
that after this resummation the infrared divergences in 
the diagrammatic technique are absent and consequently 
the external scale of the turbulence disappears from the 
diagrammatic equations. Therefore a scale-invariant so- 
lution of these equations can exist. 
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An equation for the quasi- Lagrangian velocity u which 
was introduced by Eqs. (0.1,3.2) can be derived from the 
Navier-Stokes equation; it is 

du a /dt + VpUup - u p){u a - u 0a )j + V a P 

= i/V 2 u a + f ai V Q u a =0. (3.1) 



Here v is the viscosity, P and f are qL variables related 
to the pressure P and to the pumping force f as th e qL 
velocity u is related to v, and the quantity u 0a in (3.1) 
denotes u a (t, tq). The equation (3.1) differs from the 
Navier-Stokes equation in the terms Uo subtracting the 
sweeping at the marked point r . Note that the pres- 
ence of the marked point ro in the theory means the 
loss of homogeneity: this is the price for eliminating the 
infrared divergencies from the diagrammatic expansion. 
To find correlation functions of the Eulerian velocity v 
we should express these correlation functions via correla- 
tion functions of the qL velocity u, which is not a trivial 
task. Fortunately the simultaneous correlation functions 
of the Eulerian velocity v and of the qL velocity u co- 
incide. Note that time-dependent correlation functions 
in these sets of variables are different: moreover unlike 
the qL correlation functions the time dependence of the 
correlation functions of the Eulerian velocity v is not of 
scaling type. 

The Wyld diagrammatic expansion is formulated in 
terms of propagators and vert ices determined by the non- 
linear term of the equation ( [3.l| ). The propagators are 
the Green's function G and the pair correlation function 
F of the velocities. We will treat these propagators for 
qL variables. The G-function is the linear susceptibility 
determining the average (u a ) which arises as a response 
to the nonzero average (f a )' 



G a/9 (i ) ri > r 2 ) = -«(t»a(t,ri))/*(/ / ,(0,r a )> > (3.2) 

wher e f a is the pumping force in the right-hand side of 
( |3.l|) . The F-function is the pair correlation function of 
qL velocities: 



F a p(t,r 1 ,r 2 ) = (u a (t, n)u l a(0,r 2 )) 



(3.3) 



Note that the propagators G and F in qL variables de- 
pend separately on coordinates of the points rj, r 2 which 
is the consequence of the loss of homogeneity. However 
the simultaneous correlation function F(t = 0, i"i, r 2 ) de- 
pends only on the difference ri — r 2 since it should coin- 
cide with the simultaneous correlation function of Eule- 
rian velocities. 



We expect that the propagators (3.2 3.3) possess a scal- 
ing behavior. To establish the character of this behavior 
we can utilize the KO-41 dimensional estimations. For 



the pair correlation function (3.3) we obtain from the 
estimate (|l.3[) for the characteristic velocity 



where R is the characteristic scale. For the simultane- 
ous correlation function R coincides with the separation 
between the points ri and r 2 , and for nonzero t it is de- 
termined also by the differences ri — ro and r 2 — ro, where 
ro is the marked point. The characteristic time of F can 
also be estimated. It coincides with the turnover time 
t r ~ (i? 2 /?) 1 / 3 (it is so only in qL variables). The scal- 
ing behavior of the Green's function can be established if 
one supposes that fluctuations of u does not drastically 
change the character of the response to the external force. 
Then 



G(t,ri,r 2 ) ~ R 



-3 



(3.5) 



F(i,n,r 2 ) ~ (eR) 



2/3 



(3.4) 



where R is as before the characteristic scale. The ques- 
tion arises: can such scaling behavior be obtained as a 
solution of diagrammatic equations? 

To answer this question we should first reformulate the 
diagram technique in terms of the bare vertices but with 
the dressed propagators F and G. Then one can easily 
check that the scaling behavior of F and G determined by 
the estimates (3.4,3.5) is reproduced in any order of the 
perturbation theory. But this is not sufficient to justify 
the assertion that F and G actually possess such scaling 
behavior. The reason for this was long ago recognized in 
the theory of second order phase transitions. Reformulat- 
ing the diagrammatic series for the correlation functions 
of the order parameter in terms of the bare interaction 
vertex but with the dressed pair correlation function with 
its suitable scaling exponent, one can check that this ex- 
ponent is reproduced in each order of the perturbation 
theory. Besides one immediately encounters logarithmic 
ultraviolet divergences which arise in each order of the 
perturbation expansion. After summing, the logarithmic 
corrections generate power factors which strongly renor- 
malize the naive exponents. 

Fortunately this phenomenon does not occur in the 
theory of turbulence. As was demonstrated by Belinicher 
and L'vov p3f in qL variables there are neither infrared 
nor ultraviolet diverge nces in the diagrammatic expan- 
sion for G and F, if (3.4,3.5) are assumed. The anal- 
ogous theorem can be proved for high order correlation 
functions of u and nonlinear susceptibilities: the KO- 
41 exponents are reproduced in the diagrammatic series 
and both ultraviolet and infrared divergences are absent. 
This property is the ground for the assertion that in the 
consistent theory the simultaneous correlation functions 
actually have naive KO-41 exponents. To be more pre- 
cise we should note that this property is true for the 
correlation functions of u in k-representation. In passing 
to r-representation we encounter infrared divergences re- 
lated to the fact that the average value (u 2 ) (coinciding 
with (v 2 )) is determined by the largest scale (the scale of 
pumping) and is consequently given by a formally diverg- 
ing integral. To avoid this difficulty we should consider 
the correlation functions of such objects as velocity dif- 
ferences, in expressions for which the infrared divergences 
cancel. 
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The next question is: does the absence of ultravio- 
let divergences in diagrams for correlation functions of 
u mean that all objects in the theory of turbulence can 
be characterized by KO-41 exponents? Our answer is 
"no" , since ultraviolet divergences can be observed for 
more complicated objects. Namely, we will demonstrate 
that the ultraviolet logarithms immediately arise in the 
diagrams for correlation functions of powers of the ve- 
locity gradient V a vp. The simplest example of such a 
correlation function is if ee (ri,r 2 ) = ((e(ri)e(r 2 ))} since 
the energy dissipation rate e is proportional to the sec- 
ond power of VqW/3. Those logarithms as in the theory of 
second order phase transitions lead to the renormaliza- 
tion of the scaling behavior with respect to the normal 
KO-41 one that is anomalous scaling. 



B. Ultraviolet logarithms 

Let us analyze the diagrammatic expansion for the cor- 
relation function K e£ in qL variables in terms of the 
dressed propagators F and G. The first diagram for 
K £e {v\,V2) is depicted in Fig. [|, where a w avy line cor- 
responds to the pair correlation function (|3.3| ) and an 
empty circle marks the point i"i or r 2 . These points are 
designated by 1 and 2 in Fig. |l|. Since e is proportional to 
the second power of V a vp the gradients should be taken 
of both propagators attached to an empty circle. Us- 
ing the estimate (3.4) for the function F we immediately 
obtain that the diagram depicted in Fig. [j] gives the ex- 
pression possessing the normal KO-41 behavior oc i? -8 / 3 . 
This diagram reproduces Golitsin's result |32|. No any 
divergence is produced by this diagram. To observe di- 
vergences we need to examine higher order diagrams. 



the Green's function (3.2): The wavy part corresponds 
to the variation of the velocity u and the straight part 
corresponds to the variation of the force f. The vertices 
on the diagra ms a re determined by the nonlinear term in 
the equation (3.1): there are two wavy and one straight 
line attached to each vertex. Two extra two-loop dia- 
grams should also be taken into account which can be 
obtained from the diagrams reproduced in Fig. [2] by con- 
verting left and right sides. Using the estimates ( |3.4|[3~5| ) 
we find that all of these diagrams give us expressions 
for K ee which behave as i? -8 / 3 . However there are also 
logarithmic divergences in these diagrams related to the 
loops marked by the letter "r" . The lines of the loop of 
the first diagram correspond to the product VGVGF. 
In principle there are also V's in the vert ices but the 
structure of the qL vertex determined by (3.1) is such 
that it produces the gradient of the attached propagator 
with the smallest wave vector [^3[^4), that is with the 
largest characteristic scale. This means that for loops 
with separations r (see Fig. ^|) which are smaller than 
the separation R between the points ri and r 2 , the gra- 
dients producing by the vertices are external to the loops. 
Thus the expression corresponding to the loop can be ob- 
tained by integration of VGVGF over the two times and 
coordinates of two points correspon ding to the vertices of 
the diagram. Using the estimates (3.4 J3~5) we conclude 
that this integration is dimcnsionless and produces conse- 
quently a logarithm. The upper limit in this logarithmic 
integration is R and the l ower lim it is the viscous length 
X] at which the estimates (3.4.3.5) fail, therefore the log- 
arithm produced by the right loop is \a.(Rjr\). Then the 
left loop gives i? -8 / 3 and the final contribution of the di- 
agram is oc R- & / 3 ln(R/r]). The same is correct for other 
diagrams of the same order. 




FIG. 1. The first diagram for K EE producing the normal 
KO-41 scaling. 




FIG. 2. The first diagrams for K £e producing ultraviolet 
logarithms. 



The diagrams of the next order are depicted in Fig. |2J, 
where the new combined wavy-straight lines designate 



R2 p 

p Ri 



FIG. 3. The two-loop diagram for K EE producing the sec- 
ond power of the ultraviolet logarithm. 

Let us now consider higher order diagrams for K ee . 
An example is presented in Fig. || where a three-loop 
diagram is depicted. Based on the above analysis one 
can readily recognize that this diagram gives the sec- 
ond power of the ultraviolet logarithm. The first loga- 
rithm is produced by the right loop (marked by "i?i"). 
It originates in the integration over separations R\ be- 
tween the viscous scale 77 and the characteristic separa- 
tion R 2 of the middle loop, this logarithm being \n(R 2 /rj). 
The middle loop produces an extra logarithmic integra- 
tion over dln(i? 2 /?7). The result of this integration is 
J d\n(R 2 /r])lii(R2/'n)= (l/2)\n 2 (R/rj) where R is the 
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separation between the points ri and r%. Thus we ac- 
tually find the second power of the logarithm, arising as 
a prefactor of i?~ 8 / 3 produced by the left loop. Note 
that some three-loop diagrams will not produce the sec- 
ond power of the logarithm. An example of such a dia- 
gram is given in Fig. ^. This diagram contains crossed 
lines. It prohibits the existence of the region of integra- 
tion R ^> i? 2 ^> Ri producing the second power of the 
logarithm. 




FIG. 4. The two-loop diagram for K ££ producing only the 
first power of the ultraviolet logarithm. 

Generalizing the above analysis we can conclude that 
diagrams of the n-th order will produce the normal KO- 
41 factor i? -8 / 3 with prefactors which are different pow- 
ers of the logarithm ln(i?/n) up to the n-th power. Thus 
we encounter a series over hx^R/rj), which could be an 
arbitrary function of ln(i?/n). In the next subsection we 
will argue that this function is an exponential one, which 
is a power of R/r). Such a function in the prefactor leads 
to the substitution of the normal KO-41 dimensional- 
ity by other ones: this is the mechanism producing an 
anomalous dimension. 

Let us recognize that the origin of the ultraviolet log- 
arithms (leading to anomalous dimensions) in the dia- 
grammatic series is related to the fact that instead of the 
correlation functions of the velocities we have taken the 
correlation function of the powers of the velocity gradi- 
ents. Indeed, we have seen that the logarithms are pro- 
duced by loops with additional differentiation of propa- 
gators. In our example these additional differentiation 
were related to the structure of e which is proportional 
to the second power of the velocity gradients. The loops 
without additional differentiations do not produce ultra- 
violet logarithms due to the main property of the qL 
vertex: the differentiation at the vertex is taken of the 
propagator with the largest characteristic scale. There- 
fore we should differentiate the propagators which are 
always outside a "ultraviolet loop" (that is a loop with 
a small separation). The same property of the qL vertex 
supplies the "reproduction" of extra differentiations: we 
have seen that the differentiations external to the first 
"ultraviolet loop" were internal for the next "ultraviolet 
loop" , and so on, and so forth. A repetition of such loops 
results in the accumulation of powers of the logarithms. 
It is clear that the same mechanism will work for any cor- 
relation function of objects containing velocity gradients 
(but at least two gradients). Therefore we expect that 
all such objects should possess anomalous dimensions. 



C. Anomalous Dimensions 

Let us return to the analysis of the correlation function 
K £e . In the framework of the Wyld technique a formally 
exact diagram representation for the correlation function 
can be formulated originating from the fact that in each 
diagram for K £S there exists only one cut going along 
all F-functions. This enables us to formulate the repre- 
sentation depicted in Fig. [| which is analogous to the 
representation of the imaginary parts of propagators in 
quantum field theory. Here we have classified diagrams 
for K £e in accordance with the number of F-functions in 
our marked cut (since this number runs from 1 to oo we 
obtain an infinite number of terms); the ovals designate 
objects which are sums of the diagrams at the left and 
at the right sides of the marked cut. As before empty 
circles mark the points ri and Y2 where e taken. 



cQ-x/x/x/'yQd -4- 




D -j- C 



3 + 



FIG. 5. The formally exact diagrammatic representation 
for K ee , the first terms of an infinite series. 

The first "one-bridge" term of the diagrammatic series 
depicted in Fig. ^| is determined by the block which can 
be represented as the sum of two terms depicted in Fig. [(], 
where we have introduced the new object designated by 
the oval with the inscribed line. Unlike the empty oval 
the oval with the inscribed line has one straight leg. Ac- 
tually both ovals occurring in Fig. ^ arise in the diagram- 
matic equation for the oval entering the second "two- 
bridge" term in Fig. |^. Therefore we begin our analysis 
with the equation for this oval. Diagrams for the oval 
can be classified according to the possibility to cut the 
diagram along two lines. Note that those two lines can 
correspond to two G-functions or to the functions G and 
F but not to two F-functions since the cut with all F- 
functions is unique. Designating by boxes the sums of the 
four-leg parts of the diagrams which cannot be cut along 
two lines, we come after summation to the diagrammatic 
equation presented in Fig. 0. The first term on the RHS 
here designates the bare contribution determined by the 
double V in e. Other terms are related to fluctuations. 




FIG. 6. The diagrammatic representation for the block en- 
tering the "one- bridge" contribution to K £S . 

Similar diagrammatic relations can be established also 
for the three-point object designated by the oval with 
the inscribed line and for the analogous three-point ob- 
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ject with two straight legs. The relations can be con- 
sidered as a closed system of equations for these three 
ovals. Of course the boxes entering these equations are 
not known explicitly. Nevertheless it is possible to estab- 
lish estimates for the boxes. Following the analysis given 
in (2^^Jl one can demonstrate that there are neither ul- 
traviolet nor infrared divergences in the diagrams for the 
boxes. It means that they can be estimated by the first 
contributions. The first contribution to the empty box is 
determined by the G-function with two attached vertices 
and the first contribution to the box with the inscribed 
line is determined by the _F-function with two attached 
vertices. 




33 



> + 



3 + 



3 + 




D 



FIG. 7. The diagrammatic equation for the three-point ob- 
ject entering the diagrammatic expression for K eE . 

The diagrammatic equation represented in Fig. [?] can 
be rewritten in the analytical form. For this let us intro- 
duce the designation T for the oval on the RHS of the 
diagram. Thus is a three point object and therefore T 
is a function Y a p(tx,rx, i 2 , *2, £3, 1*3), where the arrange- 
ment of the points 1, 2, 3 is shown in Fig. fj], and e is taken 
at the point 3. Then the equation represented in Fig. |?] 
takes the form 

T(£i,ri,i 2 ,r 2 ,<,r) = T (t 1 ,r 1 ,t 2 ,r 2 ,t 1 r) + 

J dt 3 d 3 r 3 dt 4 d 3 r 4 B(t 1 ,r 1 ,t 2 , r 2 , t 3 , r 3 , i 4 , r 4 ) x 



T(i 3 ,r 3 ,i 4 ,r 4 ,t, r) 



(3.6) 



Here for brevity we have omitted subscripts and an ad- 
ditional term determined by the oval with the inscribed 
line. The quantity Tq in (3.6) is the bare value of T: 



T = vS(t! - t)S{t 2 - i)V<5(r x - r)V<5(r 2 - r) . (3.7) 



The kernel B in (3.6) corresponds to the sum of boxes 
in Fig. [?] with attached lines. This kernel can be esti- 
mated as FVGVG, where the double V originates from 
t he vert ices produced by (3.1). Utilizing the estimates 
(3.4 3.5) and als o t ~ (i? 2 /e) 1//3 we conclude that the in- 
tegration in ( |3.6| ) is dimensionless. It is possible to rescale 
two other ovals to achieve di men sionless integrations in 
terms designated by dots in (O) and in two analogous 



equations for other ovals. Thus we come to homogeneous 
equations. 

Of course this homogeneity will be destroyed on scales 
of the order of r\. Therefore the structure of a solution 
will be as follows: on scales larger than ij the function T is 
a sum of power-like terms characterized by scaling expo- 
nents which could be excluded from solving the equation 
( |3.6| ) without To and with the homogeneous kernel. The 
number of these terms is infinite. This equation is lin- 
ear and consequently the coefficients of the powers cannot 
be extracted from this equation. Actually the coefficients 
can be found from the whole solution, with To taken into 
account and with account of the solution on scales of the 
order of rj. Also the set of scaling exponents does not 
depend on the details of the small-scale behavior of the 
system since it can be extracted by solving the equation 
with the kernel determined by the dynamics in the in- 
ertial interval. Returning now to K ee we conclude that 
the first two terms of the expansion presented in Fig. ^ 
produce a complicated scaling behavior of this function: 
on scales larger than rj K ee is a sum of an infinite num- 
ber of power-like terms with the universal set of scaling 
exponents but with coefficients sensitive to the dynamics 
on scales of the order of rj. All of these terms possess 
anomalous dimensions related to the scaling behavior of 
T. The same assertion can also be proven for higher or- 
der terms of t his expansion, since for all ovals equations 
of the type of (|3.6|) can be formulated leading to the same 
conclusions for the scaling behavior. 

Actually the situation is the same when we consider 
more complicated correlation functions of fields con- 
structed from velocity gradients since the gradients will 
produce logarithms and consequently anomalous dimen- 
sions. Any correlation function of such fields is as above 
a sum of terms with different anomalous dimensions. 
Therefore a question of great interest is the calculation 
of the set of anomalous dimensions. Unfortunately this 
is impossible to do analytically as we deal with a theory 
where there is no small parameter. This means for ex- 
ample that we c annot find an explicit expression for the 
kernel B in (|3.6|) ( it is possible to do for the simple model 
considered in J5l|). Therefore the set of anomalous di- 
mensions should be extracted from experiment or numer- 
ics. Nevertheless we can establish a number of selection 
rules based on symmetry reasoning. These rules allow us 
to establish a set of relations between different correla- 
tion functions of velocity gradients and to predict some 
intermediate asymptotics of the correlation functions of 
the velocity differences. These rules will be discussed in 
the next section. 



IV. FUSION RULES AND OPERATOR 
ALGEBRA IN TURBULENCE 

Here we are going to formulate for hydro dynamic tur- 
bulence a set of fusion rules for fluctuating fields as in- 
troduced by Polyakov pq, who investigated the scaling 
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behavior of the correlation functions of the order param- 
eter near a second order phase transition point. We think 
that Polyakov's idea (which is actually a kind of multi- 
pole expansion) can be successfully applied to any system 
with strong fluctuations where a scaling behavior is ob- 
served. Our treatment will be particularly based on the 
results obtained in the preceding sections but the fusion 
rules themselves can be constructed practically indepen- 
dently starting from symmetry arguments. A convenient 
language for formulating fusion rules is the so-called oper- 
ator algebra introduced by Kadanoff and Wilson [ p0| . 
In this section we will restrict ourselves to the analysis 
of simultaneous correlation functions of the velocity gra- 
dient which coincide for both Eulerian and qL velocities. 
These correlation functions, in contrast to different-time 
correlation functions of qL velocity, possess space homo- 
geneity which essentially simplifies the analysis. 



A. Scaling Dimensions of Local Fields 

Here we will treat correlation functions of local fields 
<Pj (r) constructed as different single-point products of the 
velocity derivatives. An example of a local field is the en- 
ergy dissipation rate e which is pro port ional to the sec- 
ond power of the velocity gradient 



We have seen 

in Section III that the correlation function K ee = ((ee)) 



hi) 



contains an infinite number of terms with different scaling 
behavior. The same is actually true for any correlation 
function ((<^j(r)<^j(0))) of the local fields. To proceed in 
the analysis of their scaling behavior it is worthwhile to 
extract a set of local fields A n with "clean" scaling behav- 
ior which are complicated linear combinations of single- 
point products of the velocity derivatives. The fields \&i 
and ^2 introduced by can be considered as the first 
step in this direction. Of course it is impossible to con- 
struct explicitly the final expressions for A n ; fortunately 
the expressions are not really needed for the analysis of 
the scaling behavior of correlation functions in which we 
are interested. 

Each local field A n by definition is characterized by its 
scaling dimension A„ which means that the correlation 
function (A n (R)A n (0)) cx i?~ 2A ". The question arises 
as to the scaling behavior of the correlation function 
(A n (R)A m (0)) . To establish this behavior one might re- 
call, e.g., the diagrammatic expansion for K £e presented 
in Fig. ||. The scaling behavior of the terms of this ex- 
pansion is characterized by the anomalous dimensions 
of the corresponding diagrammatic objects designated 
by the ovals. The anomalous dimension of each term 
is equal to the sum of the anomalous dimensions deter- 
mined by the pair of ovals. The same situation takes 
place for the correlation function (A n (R)A m (0)): in the 
diagrammatic language it is the pair of diagrammatic 
"fur-coats" dressing the points R and and a number of 
"bridges" made from F'-functions connecting these "fur- 
coats" . Each "fur-coat" has its own anomalous dimension 



and therefore the anomalous dimension of (A n (R)A m (0)) 
is equal to the sum of the anomalous dimensions of the 
"fur-coats" (since the F-function has no anomalous di- 
mension). Thus we come to the conclusion that the scal- 
ing index of the correlation function (A n (R)A m (0)) is 
equal to A„ + A TO . This means that 



(A n (K)A m (0)) oc R 



— A n — A. 



(4.1) 



This is the basic property of the fields A n which will be 
exploited further. 

Among the set A n of the local fields with definite scal- 
ing dimensions one can extract a subset of the so-called 
primary fields O n which give rise to all other fields A n 
by differentiation [b2| . These "field-descendants" A n are 
usually referred as secondary fields. The dimension A„ of 
any secondary field A n differs from the dimension A m of 
the corresponding primary field O m by an integer num- 
ber I: A„ = A m + I, the number I being the number of 
the differentiations needed to obtain A n from O rn . It is 
obvious that the set A m of the scaling dimensions of the 
primary fields determine also the whole set A n . 

Any local field tpj can be expanded in a series over the 
fields A n with some coefficients <Pjt n y. 



](p j{n) A n (r). 



(4.2) 



This expansion enables one to reduce the correlation 
functions of the field tpj to the correlation functions of the 
fields A n . It is convenient to order the fields A n over the 
values of their scaling dimensions: Ai < A2 < A3 . . .. 
Unfortunately it is impossible to find the values of A„ 
(to do this we s hou ld find eigenfunctions of an equation 
of the type of (3.6), the kernel of which is not known 
explicitly) but it is possible to express the scaling behav- 
ior of observable quantities in terms of A„. It is clear 
that the principal scaling behavior of correlation func- 
tions of a local field ipj is determined by the first nonzero 
term of its expansion into a series over A n . For exam- 
ple if the first terms of the expansions of (p% and if2 are 
not equal to zero the scaling behavior of the principal 
term in the correlation function (<£>i(R)y>2(0)) is related 
to the correlation function (A\Ai) since in accordance 
with the definition all other correlation functions (A m A n ) 
decrease more rapidly than (A1A1) as R — > 00. Thus 
(ipi(R)ip 2 (0)) oc R~ 2Al as R — > 00. Since the dimensions 
of the secondary fields are larger then the dimensions of 
the primary fields we first should be interested in the 
terms of the expansion of (fj over O n . 

The introduction of the fields A n enables one not only 
to predict the scaling behavior of the pair correlation 
functions but also to investigate the asymptotic behavior 
of more complicated correlation functions. Consider for 
instance a three-point correlation function of the fields ipj 
which using (4.2) can be reduced to three-point correla- 
tion functions of the fields A n . Suppose now that among 
the points there are two nearby points and the third point 
is separated from the pair by a large distance R, and we 
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are interested in the behavior of this correlation function 
as R varies. It is obvious that in this case the product of 
the fields A„(ri)A m (r 2 ) taken at the nearby points be- 
haves like a single-point object, which can be expanded 
into a series over A n (r), the point r being located near 
the points ri and r 2 . Thus we come to the relations 



^„(n)^ m (r 2 ) 



i 



C m „,i(ri -r,r 2 -r)A t (r), (4.3) 



which arc known as the operator algebra P9|j30| . Actu- 
ally it is a kind o f mu ltipolc expansion. Let us explain 
how the relation ( |4.3| ) can be found explicitly. For this 
we should expand A n (ri) and A m (r 2 ) in a series over 
ri — r and r 2 — r. Then we arrive at products of the local 
fields taken at the point r which can be reexpanded into 
a series over A n (r) 

The relation ([L3J) can be used in investigating any cor- 
relation function of the fields ipj with two nearby points 
since it means that the correlation functions of the prod- 
ucts of the left-hand and of the right-hand sides of this 
relation with any remote object coincide. For this inves- 
tigation we should first expand the prod uct tfji pj t aken 
in nearby points in accordance with (4.2), then (4.3) en- 
ables one to present the result in the form of a series 
over A n . Of course this procedure is very complicated, 
but actually we can restrict ourselves to the first terms of 
the expansion over A n which considerably simplifies the 
procedure. 



B. Local Fields in Turbulence 

In this subsection we consider the local fields (fj built 
up from the velocity derivatives V a vp from the point 
of view of symmetry. The symmetry imposes some re- 
strictions on the mutual correlation functions of various 
local fields. For example the irreducible cross correla- 
tion function of a scalar and of a pseudoscalar should be 
equal to zero. The point is that the correlation function 
should change its sign under inversion while it is impossi- 
ble to construct such a pseudoscalar correlation function 
from the separation vector R. Let us stress that differ- 
ent behavior of the fields ipj under time reversal t — ► — t 
does not impose any restriction on the type of nonzero 
correlation functions. This is because the time reversal 
symmetry in the turbulent system is violated (due to the 
presence of a nonzero energy flux) . 

Local fields with different transformation properties 
arise already in the first order in velocity gradients, 
namely the stress tensor s a p — |(V a ^ + X7pv a ) and the 
vorticity u> = V x v. The number of objects of the second 
order in V a vp is larger: there are two scalar fields s 2 and 
a; 2 , one pseudovector field s a pcop, three traceless symmet- 
ric tensor fields s ai s 1 p — (l/3)s 2 S a p, u: a ujp — (l/3)w 2 <5 Q /3 
and e a /3jU)f3S- f s + £8/3~/<-<Jf]S-y a , one irreducible pseudotensor 
of the third rank and one irreducible tensor of the fourth 
rank. It is obvious that the number of different fields 



which can be built up from V ' a vp is infinite. A field tpj 
with given transformation properties can be expanded 
into the series (4.2) over the fields A n with the same 
transformation properties: a vector is expanded into a 
series over vectors A n , an irreducible tensor of the sec- 
ond rank is expanded into a series over tensors A n _ a p 
and so on, where the coefficients in these expansions are 
scalars. 

A special consideration is needed for the correlation 
functions of the first powers of the strain tensor s a p and 
the vorticity u>. These correlation functions can be re- 
duced to gradients of the correlation function of the ve- 
locity. As we noted there are no ultraviolet divergen- 
cies related to the velocity itself. This means that both 
the strain tensor s a p and the vorticity u) have the nor- 
mal KO-41 scaling exponents 2/3. Additional restrictions 
are related to the incompressibility condition V • v = 0. 
For example, the cross correlation function of the veloc- 
ity itself with any scalar field tpj is equal to zero. To 
prove this, note that the correlation function (v(r)y>j(0)) 
is a vector which can only be directed along r. How- 
ever we know that the divergence of this vector should 
be equal to zero because of incompressibility. This is pos- 
sible only if this vector is equal to zero (the possibility 
(v(r)c^ (O)} oc r/r 3 is actually also excluded because of 
the singularity at r = 0). This means that the correlation 
functions (s a p(r)(pj(0)) and (u!(r)ipj(0)) are also zero 
since they can be obtained from (v(r)<pj (0)) by direct 
differentiation. Analogous reasons enable us to establish 
the form of the correlation function (V 7 W a (r)A n)| g(0)) 
where A n> p is a vector field with the definite scaling di- 
mension Ai jn . This correlation function should have the 
scaling exponent Ai „ + 2/3 which gives two possible ten- 
sor structures. Using now the property V • v = we find: 



(R)A n ,p(0)) 



(4.4) 



Sat 



1 - 3Ai, ra R a Rp \ Rl / 3 _ Al 



7-3Ai„ R 2 



We do not see any reason for the fields s 2 an d oj 2 to 



have nonzero first coefficients in the expansion (4.2). It 
means, for example, that the principal scaling behavior of 
all correlation functions {{s 2 s 2 }}, ((uj 2 s 2 )) and {{lo 2 lu 2 )) 
is the same: oc R~ 2Al where Ai is the scaling exponent 
of the field 0\ which is the scalar field with the smallest 
exponent. The same is true for scalar fields proportional 
to higher powers of V a vp, say (s 2 )",n > 1, (ui 2 ) n ,n > 1, 
s 3 = s a psp^s^ a , uj a s a pujp, etc. Comparing the principal 
scaling exponent 2Ai with the exponent 8/3 — 2A (where 
A is the anomalous dimension of e introduced in Section 
||) we conclude that Ai = 4/3 - A or Ai = fi/2. Be- 
sides the main terms oc 0\ in the expansion of a scalar 
field we can take into account the next term oc 2 . Such 
terms will produce corrections to the two-point correla- 
tion functions of scalar fields which scale as R~ Al ~ A2 . 

Correlation functions of vector fields or of irreducible 
tensor fields have scaling behavior different from th at o f 
scalar fields since they are expanded into a series (4.2) 
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over other fields A n , For example the principal term of 
(<Pi, a p(H) < Pj,ot0(P)) behaves oc R~ 2 - 1 , where A 2> i is the 
scaling exponent of the 2-rank tensor field Ai^ a p with the 
smallest exponent . Si nce the field s a fj is among the terms 
of the expansion (4.2) a term with a different field will be 
the leading one only if A 2i i < 2/3. In the opposite case 
the principal scaling behavior of the correlation functions 
will be determined by the s a( g-proportional term which 
means A 2 ,i = 2/3. Thus we conclude that A 2 ,i < 2/3. 
Similarly a correlation function of pseudovector fields 
scales as R~ 2Al1 with the exponent A* l < 2/3 s ince the 
field u) is among the fields A n in the expansion ( |4.2| ) for 
pseudovector fields, ft is possible to establish also a re- 
striction for the correlation functions of vector fields. The 
pair correlation functions of these fields scale as R~ 2Al 1 
with the exponent A^i < 1 + Aj since VOi is among 
the fields A n in the expansion (4.2) for vector fields. 

Note that if the system possesses conformal symmetry 
then there exists a set of strong select ion rules for the 
coefficients in the right-hand side of (4.1), established 
by Polyakov |5^]. Namely these coefficients are non-zero 
for different values A„ and A m only if these fields are 
secondary fields of the same primary field. This is the 
consequence of the "orthogonality rule" : the correlation 
functions of different primary fields O n are equal to zero if 
the system possesses conformal symmetry. This "orthog- 
onality rule" could be a basis for experimental checking: 
is fully developed hydrodynamic turbulence conformal or 
not. This question arises in connection with the recent 
work of Polyakov (53) who treated 2d turbulence in the 
framework of the conformal approach (as is known ]3l| ] 
for 2d systems conformal symmetry permits one to es- 
tablish many properties of the correlation functions, par- 
ticularly possible sets of dimensions A n ). 

One of the possible ways to check for conformal sym- 
metry in turbulence is to study experimentally the cross 
correlation fun ction ^s 2 (R)i J 9 :)Q , / 3(0))). Let us make the 
expansion ( |4.2| ) for both fields. In the case of conformal 
symmetry the nonzero contribution to the correlation 
function will produce only averages of the fields belong- 
ing to one "family" originating from some primary field. 
We expect that the main contribution will be produced 
by an average of the form (Oi(3V Q V l aOi — <5 Qi aV 2 Oi)) 
which gives 



« S 2 (R)^ a/3 (0)))oc J R-( 2 + 2A ^) 



(4.5) 



I f co nformal symmetry does not hold the exponent in 
(4.5) will not coincide with 2 + 2A 2i i. 



C. Fusion Rules for Velocity Differences 

Here we will examine the asymptotic behavior of cor- 
relation functions of the velocity differences. Consider 
the case when there are two sets of nearby points sep- 
arated by a large length R and examine the behavior 
of different correlation functions at varying R. We will 



demonstrate that in such a procedure the anomalous di- 
mensions should be revealed. Of course our treatment is 
restricted to scales belonging to the inertial subrange. 

In the preceding subsection we have considered local 
fields tfj built up from the velocity derivatives, since cor- 
relation functions of these fields have no infrared contri- 
butions. Actu ally this property enabled us to develop 
the expansion (4J2) into a series over the fields Aj with 
definite scaling dimensions. Here we wish to note that 
the behavior of the correlation functions of the velocity 
differences can be extracted from one of the correlation 
functions of the velocity derivatives using the simple re- 
lation 



1 

e (ri) - v a (r 2 ) = J drpVpVo 



(4.6) 



where the integral is taken along any curve connecting 
the point 1 to the point 2. 



The operator algebra (4.2) together with (4.6) enables 
us to reduce any product of velocity differences taken 
at nearby points to a single-point object. Consider as 
an example the se cond p ower of the velocity difference 
(vi — V2) 2 . Using (4.3,4.6) we can "fuse" this object into 
a single point r at which the expansion over all fields is 
generated. It is natural to expect that the principal term 
in this expansion is determined by 0\. 



(vi - v 2 ) 2 -> / (1) (n 2 )0 1 ((r 1 + r 2 )/2) , 



(4.7) 



where r is chosen to be equal to (ri + r 2 )/2 and r X2 
|ri — r 2 |. It means, for example, that 



«(vi -v 2 ) 2 (v 3 -V4) 2 }} cxR 



-2Ai 



(4.8) 



where ri, r 2 and r3, r 4 are pairs of nearby points sepa- 
rated by the large distance R. The r-dependence of the 
function /W (r) can also be established if one remembers 
that the general scaling behavior of the correl ation func- 
tion of velocity differences in the left-hand of (4.8) is de- 
termined by the conventional KO-41 index —4/3. Com- 
paring this index with the scaling behavior (4.£) we con- 
clude that f (1) {r) oc r Al+2 / 3 . Note that Ai+2/3 = 2-A 
where A is the anomalous dimension introduced in the 
preceding sections. 

Besides the terms of the expansion of (vi — V2) 2 over 
the scalar fields we should take into account also the 
terms of the expansion over the vector and tensor fields. 
The first term of this expansion is 

( Vl - v 2 ) 2 -> /W(r x - r a )Oi, a ((ra + r 2 )/2) , (4.9) 

where the vector fW(r) is directed along r and /i (r) oc 
r 2n/3+Ai,i ]\j t e that the terms of the expansion of (v x — 
v 2 ) over pseudoscalars or over pseudovectors are absent 
since they are forbidden by the inversion symmetry. The 
terms (Lq,4.9) give us at R 3> fi 2) 
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(((vi-v^va-v^ 2 )} 
/ (1) (ri 2 )/ (1) (r 3 4)(Oi(R)Oi(0)> 

+/P(ri 2 )/ (1) (r 3 4)(Oi,a(R)Oi(0)) 
+/ (1) (r-i 2 )4 1) (r-34)(Oi(R)Oi, 7 (0)) 

+ ^ 1) (r-12)4 1) (^34)(Ol,a(R)Ol, 7 (0)), 



(4.10) 



where R is the vector beginning at (r3 +r 4 )/2 and ending 
at (ri+r 2 )/2; r i2 = |ri -r 2 | and r 34 = r 3 — r4 1 . We see 
that there are three different contributions to the correla- 
tion function which behave oc i?~ 2Al , oc /£- A i- A i,i an d 
oc Rr 2Al1 respectively. These three contributions can be 
distinguished by their angular dependence on the angles 
between R and separation vectors ri — r 2 and r 3 — r 4 . 

The proposed scheme can easily be generalized for all 
even powers (vi — v 2 ) 2n . We again expect that the prin- 
cipal term arising as a result of "fusion" is 



(v x - v 2 ) 2 " - (ri^Or ((ri + r a )/2) 



(4.11) 



where r 12 = — r 2 |. The expression (4.11) leads to the 
conclusion that the asymptotic behavior of any correla- 
tion function (((vi — v 2 ) 2rl (v3 — V4) 2m )) is the same as 
(4.8). The scaling behavior of the function f^ n '(r) is as 
follows: f( n '(r) oc r 2 ™/ 3 + A i. Therefore we can establish 
the dependence of (((v x — v 2 ) 2n (v3 — v 4 ) 2m }} not only 
on large separations but also on small separations. 

Then we should take into consideration the terms of 
the expansion of (vi — v 2 ) 2 ™ over vector and tensor 
fields. The principal term of the expansion of (v x — v 2 ) 2 ™ 
over the vector fields is as ( |4.9[ ). This term produces 
for the correlation funct ion ( ((vi — v 2 ) 2 "(v 3 — v 4 ) 2m )) 
the same structure as ( 4.10|) . In principle one could 
take into account also the terms of the expansion of 
(vi — v 2 ) 2 ™ over high-order tensorial fields. These terms 
would be relevant if the scaling dimension A/^i of a Re- 
order tensorial field is smaller than A^i. In this case 
the angle-dependent contribution to the correlation func- 
tion (((vi — v 2 ) m (v3 — V4)™)) will be determined by the 
smallest value A ; i of Ak,i- Note that A ; i < 2/3 since 
A 2j i < 2/3. The main contributions to the correlation 
function (((vi — v 2 ) 2 " (v3 — v 4) 2m )) in this case can be 
represented like (4.10). The first contribution oc R~ 2Al 



will not depend on the angles, the second contribution 
oc sum of two terms depending on the 

angle between R and ri — r 2 or on the angle between R 
and r3 — r 4 only and the last term oc i? _2A;1 depends 
on both angles. This is again a point where a conformal 
symmetry would reveal itself : if the system possesses this 
symmetry then the contribution oc Al— A;1 is absent. 

Above we have considered the even powers of the ve- 
locity differences. Let us now analyze the correlation 
functions of odd powers. First we consider the special 
case of the first power since the difference v x — v 2 pos- 
sesses the normal KO-41 dimension. The relation (4.6) 
shows that the main term of the expansion of this differ- 
ence in a series over local fields is \J a vp. This means, for 
example, that ({v\ a — V2 a ){v3/3 — U4/3)) oc i?~ 4 / 3 . Consider 



now the correlation function ((v± a — u 2ct )(v3 — v 4 ) 2 ™). As 
we have seen the correlation function (vO„) is zero for 
any scalar field O n . Therefore only the term of the type 
(4.9) should be taken into account in the expansion of 
(v 3 



v 4 ) 2n , giving 



<(^ Q -« 2Q )(v 3 -v 4 ) 2 "}ocir 2 / 3 - 



-A.- 



(4.12) 



where A ; i as above is the smallest exponent of ten- 
sor fields entering also the correlation function (((vi — 
v 2 ) m (v 3 — v 4 )™)y. The vector structure of the corre- 



lation function ( 4.12 ) is determined both by the vec- 
tor R and the vectors ri — r 2 and 1 -3 — r 4 . Of course 



among the fields A n in the expansion (4.2) for (V3 — v 4 ) 2n 
there is a term with s a p. This means that in any case 
there is a term oc i?~ 4 / 3 in the correlation func tion 
((via ~ w 2 a)(v 3 - v 4 ) 2 "). If A ;1 < 2/3 then (gl|) de- 
creases slower with increasing R. This is again a point 
where conformal symmetry could be checked: it admits 
only the behavior oc i?~ 4 / 3 . 

Now consider a general odd power of the velocity dif- 
ference (via — w 2ct )(vi — v 2 ) 2 ". The first terms of its 
expansion into a serie s over th e fields A n has actually 
the same structure as ( 4.11 ,4~9|): 



\2n 



(Via ~ U2o)(vi - V 2 J 

/(")(r 12 )0 1 (r) + /<; ) (n 2 )Oi, /3 (i 



(4.13) 



Here f ( a n \r) oc r (2«+i)/3+A 1 ) /<W ( r ) cx r^n+D/3+A^ _ 

Analogously the terms of the expansion of (v\ a — 
w 2q )(vi — v 2 ) 2n over tensorial fields can be introduced. 
We see that the terms of the expansion of an odd power 
of the velocity difference are expressed via the same fields 
as the expansion of even powers. Therefore the behavior 
of the mutual correlation functions of the odd-odd and 
of the odd-even correlation functions at large separations 
will be the same as the behavior of the even-even correla- 
tion function. Terms with different scaling exponents can 
in principle be separated on the basis of their angular de- 
pendence which for the odd powers is more complicated 
than for the even ones. 



CONCLUSION 

We have investigated the scaling behavior of the cor- 
relation functions of the turbulent velocity in the iner- 
tial subrange of scales. We have argued that there is no 
physical reason for the correlation functions of the veloc- 
ity differences to deviate from the behavior predicted by 
Kolmogorov's dimensional estimates in the limit Re^ oo. 
This conclusion is confirmed by a rigorous treatment 
in the framework of the Wyld diagrammatic technique 
beginning with the Navier-Stokes equation: after the 
Belinicher-L'vov resummation both infrared and ultravi- 
olet divergencies are absent for the Kolmogorov solution 
determined by (|l.3| ). However correlation functions of 
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powers of the velocity derivatives possess an anomalous 
scaling behavior even in the limit Re— * oo. Their scaling 
exponents differ from those predicted by simple KO-41 
values. We present the physical picture which leads to 
anomalous exponents: they appear as a result of the tele- 
scopic multi-step eddy interaction producing anomalous 
contributions to the correlation functions. This physi- 
cal picture can also be reproduced on the diagrammatic 
level where we have extracted a series of logarithmically 
diverging diagrams whose summation leads to renormal- 
ization of the normal KO-41 dimensions. An infinite set 
of anomalous dimensions should arise. Thus the situation 
appears to be analogous to that in the theory of second 
order phase transitions. 

In order to describe the scaling behavior of correlation 
functions of local fields built up from the velocity deriva- 
tives a set of fields A n with definite scaling exponents 
A n were introduced. Any local field can be expanded 
into the series over the fields A n . The principal scaling 
behavior of the correlation functions is determined only 
by the first terms of this expansion. The symmetry clas- 
sification of the fields A n enables one to predict some re- 
lations between different correlation functions. We have 
formulated also some restrictions imposed by the incom- 
pressibility condition and propose some tests enabling 
the experimental testing of conformal symmetry of the 
turbulent correlation functions. We have demonstrated 
that anomalous scaling behavior should be revealed in 
the asymptotic behavior of correlation functions of ve- 
locity differences, and we have proposed a way to extract 
the anomalous exponents from experiment. 

Now some words concerning the applicability region 
of our theory. We considered the case of infinitely l arge 
Re where in our consideration the KO-41 scaling ( |l.3| ) 
for the velocity correlation functions is realized. How- 
ever deviat ions from the KO-41 values £n — rt/ 3 of the 
exponents ( |l . 1 1| ) of the structure functions (|l.2| ) are ob- 
served in experiment and numerics |50| . We be- 
lieve that this is related to the finite values of Re, since in 
our meaning a wide region of scales down to the viscous 
length r\ with a crossover behavior of turbulent correla- 
tion functions should occur. Actually this crossover may 
be observed in experiments. If this crossover behavior is 
characterized by a single scaling then all our conclusions 
could also be applied to this crossover region and only 
small corrections should be introduced into the values 
of the KO-41 exponents noted in Section IV. However if 
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this crossover behavior is related to a multifractal picture 
then our conclusions are not valid, and this case needs a 
separate consideration. 
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